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The Conlomb screening effect and the hnite-size effect snch as snrface tension are £g- 
nred ont in the hadron-qnark deconhnement phase transition. We stndy the mixed phase 
of the qnark droplets immersed in hadron matter. We see that the droplet phase is me¬ 
chanically nnstable if the snrface tension is strong enongh. Once the Coulomb potential is 
properly taken into account, we could effectively satisfy the condition for charge chemical 
equilibrium in the Maxwell construction. As a result, we suggest the Maxwell construction 
revives the physical meaning effectively. 

1. Introduction 

Nowadays, it has been believed that hadron matter becomes quark matter in high- 
density region, that is to say “deconhnement phase transition”. We are here interested in 
how hadron matter changes to quark matter at zero temperature. 

In the hrst-order phase transition, such as the water-vapor phase transition, we know 
the the Maxwell construction (MC) works to derive the equation of state in thermody¬ 
namic equilibrium. When we naively apply MC to this deconhnement phase transition, 
we immediately hnd one of the Gibbs conditions (GC) is not fulhlled, which are basic 
conditions in thermodynamics; 

quark hadron/_ . \ . quark _/ . hadron oQuark ohadron Tnquark hadron /i \ 

hn = hn l=hBj, /ie U he , = r , = I 

where /ig and he baryon-number and charge chemical potentials, respectively. 

Glendenning P] claimed, by the bulk calculation, that MC is not appropriate in quark- 
hadron deconhnement phase transition because of more than one chemical potential, and 
that there is assumed local charge neutrality instead of total charge neutrality; both 
hadron matter and quark matter should be charged. Consequently there appears no 
constant-pressure region like in MC and the structured mixed-phase develops in a wide 
density region. 

The bulk calculation assumes bulk quark matter and hadron matter without any hnite- 
size ehect like the surface tension and the Coulomb interaction. Therefore matter is no 
more uniform but takes geometrical structures in the mixed phase. These geometrical 
structures are called “droplet”, “rod”, “slab”, “tube” and “bubble.” 
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On the other hand, Heiselberg et al [ 2 ] pointed out that the mixed phase should appear 
in a narrow density region and it is energetically disfavored for a large value of the surface 
tension. In the previous work j2] it has been pointed out that we could satisfy the 
condition of charge chemical equilibrium even in MC case by properly introducing the 
Coulomb potential, 


fir* (2) 

where ^ which corresponds to Eq. (1) in a gauge invariant 

fashion. 

We can see that the Coulomb screening effect plays an important role and that the 
geometrical structures become mechanically unstable in the hadron-quark deconhnement 
phase transition. However, there is used the linearized approximation for the Poisson 
equation to hgure out the Coulomb screening effect analytically. 

In this study, we numerically solve the Poisson equation without any approximation. As 
mentioned above, since the Coulomb potential and its Coulomb screening effect should 
be so important to change the description of the quark-hadron mixed phase, we must 
carefully treat the Coulomb interaction. 

2. Formalism 

We use the density functional theory (DPT) to discuss the quark-hadron mixed phase 
[Sj. Consider quark droplets with a radius i?, embedded in nuclear matter. We divide 
the whole space into equivalent Wigner-Seitz cells with a radius and impose the 
total charge neutrality and chemical equilibrium in the quark and hadron phases and at 
the quark-hadron boundary. The quark phase consists of u, d, s quarks and electrons in 
chemical equilibrium. To take into account the conhnement we incorporate the bag model 
since we poorly understand the physical picture of the hadron-quark interface. We use 
the bag constant B = 120 MeV/fm^, which is equal to that in Ref. | 2 ]. As for the hadron 
phase it consists of nucleons and electrons in chemical equilibrium. We use an effective 
potential parametrized to reproduce the saturation property of nuclear matter. 

At the interface of the quark and hadron phases, we introduce a sharp boundary on 
the basis of the bag model with the surface tension a. Its value should depend on the bag 
constant and other values such as s quark mass, while we treat it as a free parameter in 
this study. 

Differentiating the thermodynamic potential D with respect to each constituent density 
pi {i = M, d, s, n,p, e) or the Coulomb potential, VcouU we obtain the equations of motion 
(EOM) for pi and for Vcoui- We numerically solve these EOM with the conditions of 
chemical equilibrium in each phase and at hadron-quark boundary; especially, we fully 
solve the Poisson equation without any approximation. Note that the Poisson equation 
becomes highly non-linear because particle charge densities are complicated functions of 
the Coulomb potential. 
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Figure 1. Density profiles pi without the Coulomb screening effect (left panel) and with 
the Coulomb screening effect (right panel). This is the case of p-Q = 1232 MeV and the 
volume fraction / = {R/R-^)^ = 0.01. Each density uniformly spreads in each phase 
and Vcoui is constant ( =0 ) in the left panel, while there is re-arrangement of the charge 
densities and Vcoui is spatially dependent in the right panel. 


3. Numerical results 

Figure 1 shows one example of the density prohles without and with the Coulomb 
screening effect. In the case without the Coulomb screening, each density is uniform in 
each phase (left panel), while it is not uniform when the Coulomb screening effect is taken 
into account (right panel). Most remarkable difference between these two hgures lies near 
the hadron-quark interface and the center of the droplet. Negatively charged particle 
densities are suppressed as they are apart from the boundary; d, s and e in the quark 
phase and e in the hadron phase. On the contrary, positively charged ones are enhanced; 
u in the quark phase and p in the hadron phase. As a result, each phase tends to be 
charge neutral by itself, which should resemble the picture given by MC. 

As has been suggested in Ref. j3], the Coulomb screening effect may induce the me¬ 
chanical instability. In Fig. 2 we demonstrate an example of the case with strong surface 
tension in our framework. We can see that the Coulomb screening effect strongly sup¬ 
presses the Coulomb energy for large R and thereby minimum of the thermodynamic 
potential disappears, which means droplet is mechanically unstable. 

We show the phase diagram for the “droplet” phase in Fig. 3, where the droplet can 
exist as an energetically favored and mechanically stable state. We also depict two curves 
denoted by “Maxwell” and “Bulk Gibbs” for comparison, which imply the one given by 
MC and the one constructed by GC under the ansatz of bulk quark matter and hadron 
matter, respectively. We can see the pb — he curve for the droplet phase appears far 
from “Bulk Gibbs”. This feature resembles the case of kaon condensation 0. While the 
droplet phase is just one of the structured mixed-phases, it could be said that even when 
there appear other geometrical structures, its he — he curve exhibits the similar behavior 
“Maxwell”. Hence, the Maxwell construction may revive the physical meaning effectively. 
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Figure 2. The thermodynamic potential 
and energy contributions as functions of 
the droplet radius R. This is the case 
of hb = 1232 MeV, / = 0.01 and a = 

60 MeV/fm^. We can see the Coulomb en¬ 
ergy is suppressed by the Coulomb screen¬ 
ing effect for large R and the minimum of 
the thermodynamic potential disappears, 
which shows the mechanical instability. 

4. Summary 

We have seen the Coulomb screening effect and the hnite-size effect in the hadron- 
quark deconhnement phase transition. We have considered here the droplet phase as one 
of the structured mixed-phases, and found that droplet is mechanically unstable when 
the surface tension is strong enough. We have found that the Coulomb screening effect as 
well as the hnite-size effect work against the structured mixed-phase and thereby restrict 
the mixed phase in a narrow density region. It could be suggested, from our work, that 
there cannot exist wide mixed phase region but only narrow region in neutron star core. 
However, we have to study other geometrical structures of “rod”, “slab”, “tube” and 
“bubble” [5] than “droplet”. After this task we can complete the phase diagram in the 
he ~ he plane. 



Figure 3. Phase diagram in the /xb — he 
plane with a = 40 MeV/fm^. In this case, 
the /iB — he curve for the droplet phase can 
appear far from “Bulk Gibbs”. Note that 
the Pb — he curve for uniform quark mat¬ 
ter denoted by “Q” is lying fairly close to 
horizontal axis because Pe ~ 0 in uniform 
quark matter. The /xb — he curve for uni¬ 
form hadron matter is denoted by “H”. 
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